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Experimental studies of Zni_ x Co x O as thin films or nanocrystals have found ferromagnetism 
and Curie temperatures above room temperature and that p- or n-type doping of Zni- x CoccO can 
change its magnetic state. Bulk Zni_ x Co x O with a low defect density and x in the range used in 
experimental thin film studies exhibits ferromagnetism only at very low temperatures. Therefore 
defects in thin film samples or nanocrystals may play an important role in promoting magnetic 
interactions between Co ions in Zni-xCoxO. The electronic structures of Co substituted for Zn 
in ZnO, Zn and O vacancies, substituted N and interstitial Zn in ZnO were calculated using the 
B3LYP hybrid density functional in a supercell. The B3LYP functional predicts a band gap of 3.34 
eV for bulk ZnO, close to the experimental value of 3.47 eV. Occupied minority spin Co 3d levels 
are at the top of the valence band and unoccupied levels lie above the conduction band minimum. 
Majority spin Co 3d levels hybridize strongly with bulk ZnO states. The neutral O vacancy and 
interstitial Zn are deep and shallow donors, respectively. The Zn vacancy is a deep acceptor and the 
acceptor level for substituted N is at mid gap. The possibility that p- or n-type dopants promote 
exchange coupling of Co ions was investigated by computing total energies of magnetic states of 
ZnO supercells containing two Co ions and an oxygen vacancy, substituted N or interstitial Zn in 
various charge states. The neutral N defect and the singly-positively charged O vacancy are the only 
defects which strongly promote ferromagnetic exchange coupling of Co ions at intermediate range. 

PACS numbers: 71.20.-b, 71.55.-i, 73.50.-h, 75.70.-i 



I. INTRODUCTION 

Following initial reports of high T c ferromagnetism 
in Zni-^Co^O 1 , there have been contradictory reports 
on the magnetic state of this system^A^. Electron 
paramagnetic resonance (EPR) measurements show that 
it is possible to produce Zni-^Co^O with x < 0.1^2^ 
in which Co is substituted for Zn rather than form- 
ing clusters. These studies also show that nanocrys- 
talline Zni-^Co^O is paramagnetic in this concentration 
range and at low temperature^; ferromagnetism was ob- 
served only at very low temperature. Another agent must 
therefore be responsible for promotion of high T c fer- 
romagnetism in Zni-^Co^O. Very recently it has been 
shown that high T c ferromagnetism can be switched re- 
versibly when Zni-^Co^O or Zni-^Mn^O nanocrystals 
are capped by O or NSiiS or when Zni-^Co^O thin films 
are exposed to Zn vapor 5 . The electric conductivity of 
ZnO changes markedly when it is annealed in vacuum 
to produce an n-type materialiiii 2 ^ while annealing in 
an oxygen atmosphere restores the insulating behavior 
of ZnO. Venkatesan et at. have shown that the magnetic 
moment per Co ion in Zni-^Co^O at room temperature 
is reduced from 2.8 \±b per Co ion to zero as the oxygen 
pressure is varied during annealling 4 . 

Here we report results of all-electron B3LYPiiii hy- 
brid density functional theory (DFT) calculations of the 
electronic states induced by defects (O vacancy, Vo, 
Zn vacancy, V ' zn-, Zn interstitial, Zn^, and N substi- 
tuted for O, No) as well as the electronic structures 
of bulk ZnO and Zni^Co^O with x = 0.028. Sev- 
eral DFT calculations on Zni-^Co^O have been reported 
recentl}ii 5 ^iSiiLi2iiS and two of thes o 16 i 18 have considered 
the role of defects in promoting ferromagnetism. One of 



the difficulties in using DFT to study defects in an in- 
sulating material with defect induced electronic states in 
mid gap is severe underestimation of the band gap and 
overestimation of exchange coupling constants 20 . In con- 
trast, the B3LYP hybrid functional predicts band gaps 
for metal oxides which are in good agreement with exper- 
imental values^ and exchange coupling constants which 
are also in good agreement with experiment. Self- 
interaction correction of DFT results in an improved pre- 
diction for the band gap of ZnO—. Zni-^Mn^O and 
Zni-^Cu^O have been studied recently using B3LYP hy- 
brid functional calculations^. 

Defects in ZnO 23 and transition metal doped 
2 n Q7 i 24 i 25 i 26 were originally studied in the 1960's and 
1970's by magnetic resonance and optical techniques. Re- 
cent magnetic resonanceSiSSi 2 ^ ^^, carrier mobility 33 
and positron annihilation spectroscopy 3 ^ and first prin- 
ciples calculations^^Si^i have given a fuller picture of 
defect induced electronic states in ZnO. Magnetic reso- 
nance techniques 8,9,38 have been applied to Zni-^Co^O 
and Zni-^Mn^O films and nanoclusters. As noted above, 
these studies 8 3*2iSi^ as well as X-ray magnetic circular 
dichroism 12,39 , which probe the local crystal structure 
about the transition metal ion, have concluded that the 
transition metal ions substitute for a Zn ion. Expitax- 
ial growth of thin films with (110) orientation have been 
achieved 12 and there is considerable anisotropy of the 
magnetism in the thin film 4 . 

The remainder of this paper is organized as follows: 
Details of calculations are described in Section II; results 
of B3LYP calculations on the defects mentioned above, 
Zni-^Co^O as well as two Co ions in the same unit cell as 
Vo, Zn^ or No are described in Section III; analysis and 
discussion of these results is given in the final Section. 
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FIG. 1: (Color online) Brillouin zone and band structure for 
bulk ZnO with the wurtzite structure. 



II. DETAILS OF CALCULATIONS 



All-electron B3LYP calculations were performed using 
the CRYSTAL code 4 ^ and a wurtzite structure (space 
group P63?7ic) with fundamental lattice constants 3.249 
and 5.206 A parallel to the a and c axes, respectively. 
Defect calculations were performed using 3x3x2 super- 
cells with dimensions 9.747 and 10.412 A parallel to the 
a and c axes. The band structure of bulk ZnO was cal- 
culated using the wurtzite primitive unit cell and 50 k 
points in the irreducible Brillouin zone (IBZ) and cal- 
culations in the 3x3x2 super cell used 13 k points in the 
IBZ. The Brillouin zone for the primitive unit cell (Fig. 
^) and the supercell have the same symmetry, hence the 
same labels are used for high symmetry points in both 
cases. All-electron Gaussian orbital basis sets were used 
for Zn 42 , O 43 , Cc 44 and N 45 . For vacancy defects a ba- 
sis set corresponding to the ion removed was included at 
the site of the vacancy ion in order to allow a proper de- 
scription of electrons localized in the vacancy. The local 
crystal structure around defects was relaxed and details 
of their structures are given in the Appendix. In most 
cases 18 ion positions, including the position of the va- 
cancy site basis set, were relaxed for each defect. These 
are ions within 3. 5 A of the defect site. Combined crystal 
structure and spin density plots were generated using the 
XCrySDen packageiML. 
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FIG. 2: (Color online) Total and atom-projected densities 
of states for bulk ZnO and ZnO with various defects from 
3x3x2 supercells. Majority spin densities of states are shown 
immediately above minority spin densities of states, (a) Total 
density of states for ZnO, (b) projected onto Co atom in Cozn, 
(c) projected onto four Zn ions surrounding Vo, (d) projected 
onto four O ions surrounding Vzn, (e) projected onto six Zn 
ions surrounding Zn;, (f) projected onto N ion in No, (g) 
projected onto N ion in No with two Co ions per unit cell, 
(h) projected onto two Co ions in No with two Co ions per 
unit cell, (i) projected onto two Co ions in Vo with two Co 
ions per unit cell. Scaling factors which have been used are 
indicated. 



III. RESULTS 

A. Bulk ZnO 

The band structure of bulk ZnO in the wurtzite phase 
is shown in Fig. Q and the total bulk density of states is 
shown in Fig. |2t- Band energies have not been shifted 
to align the valence band maximum (VBM) with any ref- 
erence energy in order to allow comparison of the band 
structure when defects are present or absent. The band 
gap of 3.34 eV obtained using the B3LYP functional is 
in good agreement with the experimental value of 3.47 
eV— and band disperison in both valence and conduc- 
tion bands is in good agreement with an earlier plane 
wave DFT study 49 and a B3LYP calculation on ZnO 22 . 
The band gap predicted for ZnO by DFT methods is 
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FIG. 3: (Color online) Band structure for Cozn in a 3x3x2 
supercell (solid line) and for bulk ZnO in a 3x3x2 supercell 
(dotted line). Left panel: majority spin. Right panel: minor- 
ity spin. 



splitting; a similar study of Zni-^Mn^G 2 ^, in which d 
levels are entirely occupied or entirely empty, shows a 
crystal field (e and £2) splitting of less than 1 eV. 

Optical spectra of Zni-^Co^O in both earlier single- 
crystal studies^2& and more recent ferromagnetic thin 
film studies^, show Co 2+ d-d transitions at mid-gap in 
ZnO. However, optical spectra in materials with highly 
localized excited states yield excitation energies for the 
N-electron state, but do not provide good estimates of 
ionization potential and electron affinity differences be- 
cause the large electron-hole attraction energy in the op- 
tically excited N-electron state reduces the excitation en- 
ergy well below the difference in ionization potential and 
electron affinity. Photoemission experiment gSSi and 
band structure calculations, which involve N+l and N-l 
electron energy levels, provide the best estimates of the 
relevant energy level positions. However, observation of 
photocurrents observed in Zni-^Co^O^ with a sub-gap 
conduction threshold may place the unoccupied d levels 
within the band gap. 



considerably underestimated, making calculations of po- 
sitions of defect levels difficult or impossible; DFT pseu- 
dopotential calculations predict band gaps of 0.23 eV 49 , 
0.91 eV 50 or 0.81 eV 35 and the gap depends strongly on 
the particular pseudopotential used^ 9 -. 



B. Substituted cobalt Coz n 

Majority and minority spin band structures are shown 
in Fig. |31and the Co atom-projected density of states is 
shown in Fig. ^jp. Occupied minority spin Co 3d states 
with e symmetry are located around 5 eV below the vac- 
uum level, while unoccupied £2 states are located around 
the vacuum level. Majority spin 3d states are fully oc- 
cupied, hybridize strongly with O 2p states and do not 
appear as a sharp feature in the density of states. By 
comparing the bulk ZnO band structure with the Co^n 
band structure in Fig. 03 it can be seen that there is a 
minimal change in the ZnO minority band structure when 
Co is introduced and that the minority spin 3d states con- 
sist of discrete levels. The majority spin conduction band 
structure of Co^n is nearly identical to that of bulk ZnO 
but there is a considerable difference between the major- 
ity spin valence band structure for ZnO and Co^ n . Re- 
sults from resonant photoemission experiments 3 ^^ show 
a well defined peak at the top of the valence band as well 
as a broader density of states due to Co extending down 
to 8 eV below the Fermi level; the photoemission peak 
near the Fermi level may therefore arise predominantly 
from emission from 3d minority spin states. The large 
energy splitting of occupied e and unoccupied t minor- 
ity spin d electrons found here is caused by differences 
in exchange energies of these levels and not crystal field 



C. Oxygen vacancy Vo 

B3LYP calculations were performed for one O vacancy 
in a 3x3x2 supercell. The deep defect level associated 
with the two 'dangling electrons' in a neutral O vacancy is 
shown in Fig. 21 and the densities of states projected onto 
the four Zn ions surrounding the vacancy and the vacancy 
site itself is shown in Fig. There is no spin polarisa- 
tion for the neutral vacancy and only a minor change in 
the ZnO valence band structure. The position of the Vo 
and Vq levels have been determined to lie 2.1 and 2.3 eV 
below the conduction band minimum (CBM) by means of 
deep level nonlinear spectroscopy 52 . The B3LYP calcula- 
tion on the neutral unit cell places the Vo e(0|+) level 3.0 
eV below the CBM; when the calculation was repeated 
with the unit cell singly-positively charged, the degen- 
erate Vq vacancy level splits into occupied e(+|2+) and 
unoccupied levels, with the occupied and unoccupied lev- 
els 3.1 and 1.6 eV below the CBM, respectively, at the T 
point. These values are somewhat deeper than those ob- 
tained from nonlinear spectroscopy, however, inspection 
of Fig. El shows that these levels have their lowest values 
at the T point. (The CRYSTAL program allows charged 
periodic unit cells to be treated by introducing an effec- 
tive, compensating uniform background charge). There 
is considerable dispersion of the defect level because tails 
of wave functions associated with the defect induced level 
extend beyond the supercell boundaries. Previous DFT 
calculations^^ on this defect have shown that the de- 
fect level wave function has large amplitude on the shell 
of O ions outside the shell of Zn ions closest to the va- 
cancy site, which may explain why this level lies so deep. 
Supercell calculations on cells larger than the 3x3x2 cell 
used are not currently feasible. 



FIG. 4: (Color online) Band structure for Vo in a 3x3x2 
supercell (solid line) and for bulk ZnO in a 3x3x2 supercell 
(dotted line). 



D. Zn vacancy Vzn 

Zn vacancies are believed to be the dominant deep ac- 
ceptor states in ZnO 34,54 . All of the charge states of 
defects discussed so far, with the exception of the singly- 
positively charged Vo state, have closed shell ground 
states and no unpaired electrons. The neutral Zn va- 
cancy contains two 'dangling holes' so that the minority 
spin band structure consists of two singly-occupied bands 
and two empty bands (Fig. [5J. These mainly consist 
of four O 2p orbitals on the O ions which surround the 
vacancy site and they form non-degenerate and quasi- 
triply-degenerate states with a and t point symmetries. 
In the neutral Vzn defect the a state lies just above the 
VBM and the t state lies above that. The degeneracy 
of the t state is broken to allow occupation of one of 
these orbitals. Breaking of the t degeneracy ought to 
lead to a Jahn- Teller distortion of the tetrahedra of O 
ions surrounding the vacancy, but there was little sym- 
metry breaking when two shells of atoms surrounding 
the vacancy were relaxed. The two electrons in singly- 
occupied orbitals exist either as a spin singlet (S=0) or 
triplet (S=l). The wave function form used consists of 
single determinants for up and down spins and can only 
describe the spin triplet state of the vacancy. It is there- 
fore not possible to predict whether the neutral Vzn 
vacancy will have a non-magnetic, singlet or magnetic, 
triplet ground state on the basis of these calculations. 



E. Zn interstitial Zn^ 

Shallow donor states estimated to lie 30 meV below 
the ZnO CBM have been observed in ZnO subjected to 
electron irradiation and have been assigned to Zn^ or Zn^ 



FIG. 5: (Color online) Band structure for Vzn in a 3x3x2 
supercell (solid line) and for bulk ZnO in a 3x3x2 supercell 
(dotted line). Left panel: majority spin. Right panel: minor- 
ity spin. 



complexes 33 . The energetically favorable site for a single 
Zn interstitial in ZnO is the octahedral site^S. The band 
structure for a 3x3x2 ZnO supercell with a single Zn in- 
terstitial at the energy minimzed octahedral site is shown 
in Fig. El A new defect level is formed after hybridization 
of the As level of the Zn interstitial and the lowest conduc- 
tion band (which is predominantly of Zn As character). 
The density of states for the Zn interstitial and the six Zn 
ions nearest to it is shown in Fig. [2^ The atom-projected 
density of states plot shows that the defect level lies just 
below the CBM of bulk ZnO (shown in Fig. |S|as a dotted 
line) The defect level disperses by about leV throughout 
the Brillouin zone and the conduction band is consider- 
ably modified at low energies. In a larger cell with much 
reduced defect interactions, this strongly dispersive level 
would likely consist of a sharp level within a few tens of 
meV of the CBM, which would confirm the assignment 
of the state observed in experiment 33 . 

Zn^ in a positive ion state has been proposed as a pos- 
sible mediator of ferromagnetism in Zni-^Co^O 5 . The 
spin density for Zn^ in a 3x3x2 supercell with one electron 
removed per unit cell is shown in Fig. {7\ The spin den- 
sity plot shows that the spin density for the positive ion 
state is nearly completely localized on the interstitial ions 
plus the 12 Zn and O ions surrounding it; this was also 
found to be the case for the charge density for the neu- 
tral interstitial associated with the defect level just below 
the CBM. The magnetic moment on the interstitial Zn+ 
ion determined by a Mulliken population analysis is 0.37 
Hb and the moments on the 6 Zn and 6 O ions immedi- 
ately surrounding the vacancy are 0.049 (Zn x3), 0.027 
(Zn x3), 0.018 (O x3) and 0.013 (O x3) /is, giving a net 
moment of 0.69 /ig on the 13 ion Zn^ complex and 0.31 
Hb on the remaining ions in the supercell. 
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FIG. 6: (Color online) Band structure for Zm in a 3x3x2 
supercell (solid line) and for bulk ZnO in a 3x3x2 supercell 
(dotted line). 




FIG. 7: Spin density for Zn^ in a 3x3x2 supercell with unit 
positive charge. Upper panel: spin density in a plane con- 
taining interstitial Zn ions. Lower panel: spin density in a 
plane perpendicular to the c axis containing Zn ions closest 
to interstitial ions. O ions are indicated by small grey spheres 
and Zn ions by large grey spheres. 



FIG. 8: (Color online) Band structure for No in a 3x3x2 
supercell (solid line) and for bulk ZnO in a 3x3x2 supercell 
(dotted line). Left panel: majority spin. Right panel: mi- 
nority spin. The position of the Fermi level is shown as a 
horizontal dotted line in the minority spin band structure. 



F. Substituted nitrogen No 

Recently it has been shown that N doping can be 
used to produce p-type ZnO 55 and there have been sev- 
eral theoretical studies of co-doped N in Zni-^Co^O and 
Zni-^Mn^Oi^ffi. Nitrogen forms a substitutional defect 
in ZnQ2£ where it leads to introduction of a hole. Sym- 
metry constraints allow the hole to be mainly associated 
with a nearly degenerate pair of N 2p orbitals lying in the 
ab plane of the wurtzite structure or in a non-degenerate 
2p orbital parallel to the c axis. The former orbital oc- 
cupation was found in the geometry optimized structure 
for No and the semi-metallic band structure is shown in 
Fig. [51 where the partially occupied, weakly dispersive 
and fully occupied, non-degenerate N 2p states lie above 
the ZnO occupied states. 

A plot of the spin density is shown in Fig. [§1 The spin 
density from the N ion delocalizes onto nearby O ions in 
the same plane (perpendicular to the c axis) and nearly 
all spin density is confined to this plane. The magnetic 
moment on each N ion is 0.66 hb and on the O ions with 
the largest moments, it is 0.03 \±b- 



G. Substituted cobalt with V™ + , No or Zn™ + 

The possibility that neutral or positively charged de- 
fects can mediate ferromagnetic interactions in ZnO at 
intermediate range (i.e. just beyond near neighbor dis- 
tances) with a strength that could result in Curie temper- 
atures above room temperature was investigated using 
total energy calculations. A defect which mediates ferro- 
magnetism will most likely be magnetic itself as induced 
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FIG. 9: Spin density for No in a 3x3x2 supercell in a plane 
perpendicular to the c axis containing only N and O ions. O 
and N ions are indicated by small grey and black spheres, 
respectively, and Zn ions by large grey spheres. 



magnetic polarisation of the defect is unlikely to result in 
strong mediation of ferromagnetism. Total energy calcu- 
lations were performed with two Co ions and one defect 
per 3x3x2 unit cell and the magnetic moments of the Co 
ions parallel or anti-parallel. The minimum distance be- 
tween Co ions in the lattice was 7.66 A and the Co-Vo 
and Co-No minimum distances were 4.94 and 5.00 A. 
The same Zn sites were substituted for Co in Co/Vo, 
Co/No and Co/Zn^ calculations. The same O sites were 
substituted by Vo or N allowing a direct comparison of 
the efficacy of magnetic coupling of these defects in var- 
ious charge states. In Co/Zn^ calculations the minimum 
Co-Zn^ distances were 5.10 and 6.68 A. The total energy 
difference between the two magnetic states for two Co 
ions in the absence of a vacancy is less than 1 meV. When 
a Vo defect was introduced, the total energy differences 
for the two magnetic states were less than 1 meV for Vo 
or Vq + vacancies and 36 meV for a Vq vacancy (with 
the ferromagnetic state of the Co ions lower in energy). 
The spin density for two Co ions and a Vq vacancy with 
the Co majority spins parallel to each other is shown in 
Fig. [HII and the density of states projected onto the two 
Co ions is shown in Fig. When a Zn^ interstitial was 
introduced, the total energy difference was less than 1 
meV for either the neutral or singly-positively charged 
states of the defect. When a No defect was introduced, 
the ferromagnetic configuration of the Co ions was 27 
meV lower in energy than the anti-ferromagnetic config- 
uration. 



IV. DISCUSSION 

A. Single defects 

Recent positron annihilation experiments 32 found that 
the Vo defect in ZnO subjected to irradiation by high 
energy electrons is neutral at low temperature and that 
these vacancies have a deep donor character with an ion- 
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FIG. 10: Spin density for two Co ions and a Vj defect in a 
3x3x2 supercell. O ions are indicated by small grey spheres 
and Zn and Co ions are indicated by large grey and black 
spheres, respectively. Spin density on the defect appears as a 
roughly spherical density while spin density on Co ions has a 
mixed d and tetrahedral lobe character. 



ization level 100 meV below the CBM. This defect was 
found to be ten times more abundant than the Vzn va- 
cancy. Thermodynamic arguments ^§££L£L which take en- 
ergy levels from DFT calculations corrected for the band 
gap error, have shown that Vo is the most thermody- 
namically stable defect under Zn rich conditions. The 
positions of the e(0|+) level of Vo are predicted to be 
2.32i or 2.4 eV 36 below the CBM and are therefore in 
good agreement with the value of 2.1 eV determined from 
nonlinear spectroscop;^ 2 -. This level may also be involved 
in the 600 nm photoluminescence (PL) observed in na- 
tive ZnO by optical detection of electron paramagnetic 
resonance (ODEPR) 31 ; one model for this PL places the 
Vo e(0|+) level 2.5 eV below the CBM^i. Using thermo- 
dynamic arguments, Kohan and coworkers^ - and Zhang 
and coworkers^ find that the Vo defect is a 'negative 
LP defect and that the e(+|2+) level lies 0.6 eV below 
the CBM 3 ^; in this case the only charge states of the Vo 
defect which would be observed in equilibrium are the 
neutral or doubly-positively charged states and it would 
be unlikely that isolated Vo defects would be the main 
mediators of ferromagnetism in Zni-^Co^O. 

Positron annihilation studies of electron-irradiated 
ZnO 32,34 have shown that Vzn vacancies exist in a neg- 
atively charged state and have given an estimate for the 
e(2-|-) level at 2.3 eV below the CBM. Estimates of the 
position of this level from DFT calculations 50 place this 
level at 2.6 eV below the CBM. We have not performed 
calculations on the negatively charged Vzn vacancy, how- 
ever, the single particle acceptor levels (Fig. EJ) lie 2.3 eV 
below the CBM, in agreement with both the positron 
annihilation and DFT results. 

The Zn^ interstitial is known to be a shallow donor and 
estimatesSi^ 3 ^ 2 ^* 5 ^ of the position of the donor level 
range between 30 meV 33 and 200 meV— below the CBM. 
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FIG. 11: Energy level diagram for two Co ions and an oxy- 
gen vacancy in the ferromagnetic ground state. The exchange 
integrals which are important for ferromagnetic exchange cou- 
pling are shown as J symbols. 

DFT studies^Si where corrections have been made for 
the band gap error, find that this level lies above the 
CBM so that Zn^ is spontaneously doubly ionized. The 
B3LYP calculations presented here show that the Zn As 
level in Zn^ hybridizes strongly with the ZnO conduc- 
tion band and forms a strongly dispersive band which 
lies just below the CBM. A calculation on a larger super- 
cell (say 5x5x3) would allow a more reliable prediction of 
the charge state and donor level position of this defect. 

B3LYP calculations for No presented in Fig. |S|show 
a semi-metallic electronic structure for this defect with 
the Fermi level 2.2 eV below the CBM. An earlier DFT 
study of this system found the acceptor level for No to 
be 0.4 eV above the VBM— . 



B. Exchange coupling mechanism 

Substantial ferromagnetic exchange coupling has been 
found only in two instances in this work and both involve 
a defect with a magnetic moment and defect spin densi- 
ties which overlap with Co ion spin densities. The ground 
state configuration is one where majority Co t<i spins are 
parallel and minority e spins are parallel to each other 
and the oxygen vacancy spin (Fig. ITT]) , so that there 
are exchange couplings between these three spins which 
lead to an overall ferromagnetic ground state. Direct 
exchange coupling of Co ions at this separation is less 
than 1 meV, as demonstrated by total energy calcula- 
tions when no oxygen vacancy is present. When the spin 
on the defect is removed by a double ionization of the de- 
fect the exchange coupling is reduced from 36 to 1 meV, 
demonstrating the importance of the spin of the defect in 
promoting exchange coupling. No substantial exchange 
coupling was found for the positively charged Zn+ defect 
which has spin half also. This may be because the spin 
associated with the interstitial Zn ion does not overlap 
the Co ion e wave function to an appreciable extent in 
the configuration studied. 

The exchange coupling mechanism outlined here is es- 
sentially the same as the impurity band model of Coey et 
al 60 in which there is strong spatial and energetic over- 
lap between an impurity band associated with a defect 
and a filled d band. A double exchange mechanism of 
ferromagnetic coupling in Zni-^Co^O and Zni-^Mn^O 



has been outlined, for example by Petit et al 18 or Kittil- 
stved et a/. 6,61 , in which certain charge fluctuations are 
only permitted for the ferromagnetic state. These charge 
fluctuations generally lower the total energy and changes 
to the energy and wave function can be estimated using 
perturbation theory. In the approach used here, energy 
differences and changes to the wave function are obtained 
in the self consistent field calculation. Mulliken spin and 
charge populations for Co ions and for the O vacancy 
and four neighboring Zn ions in both ferromagnetic and 
antiferromagnetic configurations show that when there 
is strong ferromagnetic coupling and Co spins are paral- 
lel, the spin localized on the vacancy decreases by around 
0.1 fiB compared to when the spins are antiparallel. This 
shows that there is greater derealization of the vacancy 
spin (and consequently greater overlap with the region 
where Co e levels are located) when Co spins are paral- 
lel 

Experiment and DFT calculation^ 36 ! 50 indicate that 
Vo is the most abundant defect in ZnO under zinc rich 
conditions and that Vzn has an abundance an order of 
magnitude less than that of Vo 32 - Experiment 4 shows 
that magnetism of Zni-^Co^O in n-type ZnO is de- 
stroyed by annealing in oxygen. The calculations pre- 
sented here show that Vo must be in a singly-positively 
charged state in order for it to act as an effective promo- 
tor of ferromagnetism. Thermodynamic arguments based 
on DFT calculations^ 5 ^ have been made to suggest that 
Vo is a 'negative LP defect in which a pair of Vq defects 
will spontaneously transform into Vo + Vq~ . Obviously 
Vq must be stable if indeed a magnetic defect is neces- 
sary for stabilization of the ferromegnetic ground state 
in Zni-^Co^O. In n-type Zni-^Co^O with an excess of 
Zn we therefore expect to find both Vo and Vzn defects 
with more abundance of Vo- If some of these transfer 
a single electron to Vzn or another acceptor defect then 
the magnetic defects necessary for strong ferromagnetic 
coupling of Co ions may exist. 



C. Summary 

The electronic structures of wurtzite ZnO with Co sub- 
stituted for Zn, oxygen and zinc vacancies, a zinc in- 
terstitial and nitrogen substituted for oxygen have been 
computed using the B3LYP hybrid density functional in 
3x3x2 supercells. Total energy calculations have been 
performed for two Co ions in a 3x3x2 supercell with an 
oxygen vacancy, zinc interstitial and substituted nitro- 
gen. Neutral and positively charged states of the oxy- 
gen vacancy and zinc interstitial were studied. Only 
the singly-positively charged oxygen vacancy and neutral 
substituted N defects resulted in a substantial exchange 
coupling of two Co ions at intermediate range (> 7 A). 
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APPENDIX A: STRUCTURE OPTIMIZATION 

Minimization of the total energy of the Co zn defect, 
by varying positions of ions which lie within 3.5 A of the 
Co 2+ ion, results in changes in Co-0 bond lengths to 1.98 
A parallel to the c axis and 1.97A in the ab plane. This 
corresponds to an expansion of the Zn-0 bond parallel 
to the c axis by 0.03 A and a slight contraction of bonds 
nearly in the ab plane by 0.01 A. 

When positions of ions within 3.5 A of the Vo va- 
cancy site were varied to minimize the total energy, the 
Zn ions immediately neighboring the vacancy site relaxed 
inwards by 0.25 A for the Zn-0 bond parallel to the c 
axis and by 0.14 A for Zn-0 bonds nearly parallel to the 
ab plane so that the nearest Zn-Zn distances for these 



ions were reduced from 3.21 and 3.25 A to 2.89 and 2.97 
A, respectively. These relaxed bond distances are several 
percent shorter than those obtained previously in a DFT 
pseudopotential calculation which used a plane wave ba- 
sis set 50 but are close to the distances found for the neu- 
tral oxygen vacancy in ZnO in a more recent plane wave 
DFT studyS. 

O ions around the Zn vacancy relax outwards by 0.14 
A for the Zn-0 bond parallel to the c axis and by 0.16 
A for the remaining three O ions nearest to the vacancy 
site. 

The 12 Zn and O ions closest to the interstitial Zn, as 
well as the interstitial itself, were relaxed for the neutral 
Zn^ vacancy calculations described above. The optimized 
Zn^-Zn distances were 2.36 and 2.47 A and the Zn;-0 
distances were 2.06 and 2.68 A. 

Ionic radii of O and N ions are similar and so only 
a small distortion of the ZnO lattice is found when the 
structure of the No defect is allowed to relax. The Zn-N 
bond lengths parallel to the c axis and nearly parallel to 
the ab plane, are 1.94(6) A and 1.97(6) A, respectively 
which compares to Zn-0 bond lengths of 1.952 and 1.985 
A in the bulk ZnO. 
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